SIMPLE EXPLICIT FORMULA FOR COUNTING LATTICE 
!>' POINTS OF POLYHEDRA 
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(NJ . JEAN B. LASSERRE AND EDUARDO S. ZERON 

rS^ , Abstract. Given z e C n and A e Z mxn , we consider the problem of 

evaluating the counting function h(y; z) := z x \ x&T\ Ax=y, a;>0}. 
We provide an explicit expression for h(y; z) as well as an algorithm 
with possibly numerous but very simple calculations. In addition, wc 
exhibit finitely many fixed convex cones of R" explicitly and exclusively 
defined by A such that for any y g Z m , the sum h(y; z) can be obtained 
by a simple formula involving the evaluation of z x over the integral 
points of those cones only. At last, we also provide an alternative (and 
pH , different) formula from a decomposition of the generating function into 

simpler rational fractions, easy to invert. 
Keywords: Computational geometry; lattice polytopes. 
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1. INTRODUCTION 

Consider the (not necessarily compact) polyhedron 

g! (1.1) n(y) := {x G R n | Ax = y; x > 0}, 

t"^ ■ with y G Z m and A G Z mx ™ of maximal rank for n > m; besides, given 

z G C™, let h : Z m — > C be the counting function 



(1.2) y^h{y;z):= ^ 

x£Si(y)r\Z 



(where z x stands for Ilfc^fc 21 * 1 )- ^ ne complex vector z G C" may be chosen 
close enough to zero in order to ensure that h(y; z) is well defined even when 
rS \ Q(y) is not compact. If Q(y) is compact, then y i— > h(y; z) provides us with 

the exact number of points in the set il(y) HZ" by either evaluating h(y, 1), 
or even rounding h(y; z) up to the nearest integer when all the entries of z 
are close enough to one. 

Computation of h has attracted a lot of attention in recent years, from 
both theoretical and practical computation viewpoints. Barvinok and Pom- 
mersheim [3], Brion and Vergne [8], have provided nice exact (theoretical) 
formulas for h(y; z); see also Szenes and Vergne [15]. For instance, Barvinok 
considers z *— > h(y; z) as the generating function (evaluated at z := e c G C n ) 
of the indicator function x i— > /n(y)nz n ( x ) of the set fi(y) H Z n and provides 
a decomposition into a sum of simpler generating functions associated with 
supporting cones (themselves having a signed decomposition into unimodu- 
lar cones) . We call this a primal approach because y is fixed, and one works 

l 
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in the primal space M. n in which Q(y) is defined. Remarkably, Barvinok's 
counting algorithm which is implemented in the software LattE (see De 
Loera et al. |10j ) runs in time polynomial in the problem size when the di- 
mension n is fixed. The software developed by Verdoolaege [17] extends the 
LattE software to handle parametric polytopes. On the other hand, Brion 
and Vergne [5] consider the generating function H : C m — > C of y i— > h(y; z), 
that is, 

n 

(1.3) w ~ H( W ) : = %5*K = II I _ —A k • 

y£Z m k=l 

They provide a generalized residue formula, and so obtain h(y; z) in closed 
form by inversion. We call this latter approach dual because z is fixed, and 
one works in the space C m of variables w associated with the m constraints 
Ax = y. 

As a result of both primal and dual approaches, h(y; z) is finally expressed 
as a weighted sum over the vertices of Similarly, Beck [5], and Beck, 

Diaz and Robins [6] provided a complete analysis based on residue techniques 
for the case of a tetrahedron (m = 1). Despite its theoretical interest, 
Brion and Vergne's formula is not directly tractable because it contains many 
products with complex coefficients (roots of unity) which makes the formula 
difficult to evaluate numerically. However, in some cases, this formula can 
be exploited to yield an efficient algorithm as e.g. in [2] for flow polytopes, 
in [7] for transportation polytopes, and more generally when the matrix A 
is totally unimodular as in [9]. Finally, in [I2j[l3], we have provided two 
algorithms based on Cauchy residue techniques to invert H in (|1.3|) . and an 
alternative algebraic technique based on partial fraction expansion of H. A 
nice feature of the latter technique of [13] is to avoid computing residues. 

Contribution: Our contribution is twofold as it is concerned with both 
primal and dual approaches. On the primal side, we provide an explicit 
expression of h(y; z) and an algorithm which involves only elementary op- 
erations. It uses Brion's identity along with an explicit description of the 
supporting cones at the vertices of f2(y). It also has a simple equivalent 
formulation as a (finite) group problem. Finally, we exhibit finitely many 
fixed convex cones of M. n , explicitly and exclusively defined from A, such 
that for any y S Z m , the sum h(y; z) is obtained by a simple formula which 
evaluates Yl z x over the integral points of those cones only. 

On the dual side, we analyze the counting function h, via its generating 
function H in (jl.3p . Inverting H is difficult in general, except if an appro- 
priate expansion of H into simple fractions is available, as in e.g. |13j . In 
their landmark paper [8], Brion and Vergne provided a generalized residue 
formula which yields the generic expansion 

(1.4) gw-EEU ..^ 
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Here, a G J a whenever A a is invertible, q is the smallest common multiple 
of all | det v4 CT | ^ 0, p q = e 2n ^ q is the (/-root of unity, 5k <a G {0, 1/g}, and 
Q g ,a G C. The finite group G CT has q m elements. The coefficients Q S)CT are 



difficult to evaluate. Our contribution is to expand H in (jl.3p in the form 



(1.5) H(w 



E 



1 



ZjW A i 



1 



E 



iW W ;- 4T ?[ cr > Vl 



where: Z 
(1.6) 



i? 2 (cr; z) 

, - {0, 1, . . . , fj, a — 1}, /v = I det A^l, each 7j[a, G Z n and: 
R 2 (a;z) := [l 



s-rTrn — m 



Identity (|1.5p is a nontrivial simplification of the residue formula (jl.4[) be- 
cause the 7y[<r, tt^J's are given explicitly. And so the coefficients of the rational 
fraction (II. 5p in w; are very simple to evaluate with no root of unity involved 
(it can also be done symbolically); however this task can be tedious as for 
each a G J a one has | det A a \ n ~ m terms rj[a, u^] to determine. But once 
determined, (|1.5p is easy to invert and provides h(y; z) for any y G Z m . 



2. BRION'S DECOMPOSITION 

2.1. Notation and definitions. The notation C, M and Z stand for the 
usual sets of complex, real and integer numbers, respectively. Moreover, 
the set of natural numbers {0, 1, 2, . . .} is denoted by N, and for every nat- 
ural number fj, G N, the finite set {0, 1, . . . , /i — 1} of cardinality fj, is de- 
noted by Z„. The notation B' stands for the transpose of a matrix (or 
vector) B € M sx *; and the kth column of the matrix B is denoted by 
-Bfc := (Bx ; k, . . . , B Si k)'. When y = 0, the cone 11(0) in (jl.ip is convex, 
and its dual cone is given by, 

(2.1) O(0)* := {b G R n | ft'x > for every x G 0(0)}. 

Notice that 0(0)* = M n if O(0) = {0}, which is the case if O(y) is compact. 

Definition 1. Let A G Z mx ™ be of maximal rank. An ordered set a = 
{eri, . . . , <r m } of natural numbers is said to be a basis if it has cardinality 
\a\ = m, the sequence of inequalities 1 < o\ < 02 < • • • < cr m < n holds, 
and the square [m x m] submatrix : 

(2.2) A a := [A ai \A a , 2 1 • • • \A a J is invertible. 
We denote the set of all bases a by Ja- 

Definition 2. Given a maximal rank matrix A G Z mxn , and any basis 
cr G Ja, the complementary matrices A a G Z mxn and ^ G Z mx ( n ~ m ) 
stand for L4fc]fc G(T and L4fc]fc^o-, respectively. Similarly, given z G C n , the 
complementary vectors z a G C m and z<f G C n ~ m stand for (zfc)fceo- and 
tzk)k$<r, respectively. 
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For each basis a G I a with associated matrix A a G Z' mxm , introduce the 
indicator function 5 a : Z m — ► N defined by : 



(2.3) y h- 5 a (y) :-- 



1 if A^y el 
otherwise. 



Notice that 5 a is a multi-periodic function with periods A a and \i a 
| det A a \, meaning that 5 a (y + A a q) = 5 a (y + fj, a q) = 5 a (y) for all y,q G 2 
Finally, given a triplet (z, x,u) G C n xZ"x R n , introduce the notation 



(2.4) ||*|| 
ln(z) 



max{[zi|, |z 2 |, • • • , \z n \}, 
(ln(zi),ln(z 2 ), . . . ,ln(z n )). 



Notice that z x = z%° zJ ', for all bases a G I a and all z G C n , x G 



2.2. Brion's decomposition. Let Q(y) be the convex polyhedron in (jl.ip 
with y G Z m , ^4 G Z mxn being of maximal rank, and let h : 7L m — > C be the 
counting function in (|1.2|) . with < 1. 

Obviously z) = whenever the equation Ax = y has no solution 
x G N n . The main idea is to decompose the function h following Brion's 
ideas. Given any convex rational polyhedron P C M n , let [P] : R ra — > {0, 1} 
be its characteristic function, and f[P] : C — > C its associated rational 
function, such that 

(2.5) z ^ f[P,z] := Yl 

xePnz n 

holds whenever the sum converges absolutely. For every vertex V of P, 
define Co(P, V) C W 1 to be the supporting cone of P at V . Then, Brion's 
formula yields the decomposition : 

(2.6) [p] = Yl i c °( p ,n, 

vertices V 

modulo the group generated by the characteristic functions of convex poly- 
hedra which contain affine lines. And so, 

(2-7) f[P,z}= f[Co(P,V),z]. 

vertices V 

The above summation is formal because in general there is no z G C n for 
which the series 

^{/|xePnZ n } and ^{z x | x G Co(P, V) (~1 Z n } 

converge absolutely for all vertices V. The notation ]TJ E stands for the sum 
of all elements of a countable set E C C. It is a complex number whenever 
the resulting series converges absolutely; otherwise it stands for a formal 
series. 
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Example: Let P := [0, 1] C M so that Co(P, {0}) = [0, +oo) and Co(P, {1}) = 
(— oo, 1]. Simple enumeration yields f[P, z] = z° + z = 1 + z, but one also 
has: 

f[P,z] = /[Co(P,{0}),z]+/[(P,{l}),z] = l/(l- z ) + z 2 /(z-l) = 1 + z. 

3. Computing h(y;z): A primal approach 

Let C(Xa) := {Ac | x € N n } C M m be the cone generated by the columns 
of A, and for any basis a G Ja, let C(er) C M m be the cone generated by the 
columns A^ with k G a. As A has maximal rank, C(Ja) is the union of all 
C(cr), er G Ja- With any y G C(Ja) associate the intersection of all cones 
C(o") that contain y. This defines a subdivision of C(Ja) into polyhedral 
cones. The interiors of the maximal subdivisions are called chambers. In 
each chamber 7, the polyhedron Q(y) is simple, i.e. A~ l y > for all u G Ja 
such that vl^y — 0. 

For any chamber 7, define, 

(3.1) B(X4,7) := W e Ja I 7 C C(<r)}. 

The intersection of all C(er) with <r G B(Ja,i) is the closure 7 of 7. 

Back to our original problem, and setting P := fi(y), the rational function 
/[-P, 2] is equal to h{y;z) in (|1.2I) whenever ||z|| < 1. We next provide an 
explicit description of the rational function f[Co(P,V),z] for every vertex 
Vol P. 

Let S a be the function defined in (|2.3|) ; and let Z M(T := {0, 1, . . . , — 1} 
with //(j := I det A a \. A vector V G M n is a vertex of P = Cl(y) if and only 
if there exists a basis a G J a such that : 

(3.2) ^ = A~ x y > and = 0, 

where V a and are given in Definition [2j Moreover, the supporting cone 
of P at the vertex V is described by : 

(3.3) Co(Q(y),V) := {x£« n \Ax = y; x k > if V k = 0} . 
Let us now define the larger set 

(3.4) C(n{y), a) := {x G M n | + A/zy = j/; ay > 0}, 

so that Co(fi(y),V) is a subcone of C(f2(y),o") for all bases a G J a and 
vertex 1/ of which satisfy Ty = (recall (|3.2p ). Besides, when ly = 
and y G 7 for some chamber 7, then C(fi(?/),cr) and Co(fi(y),V) are 
identical because fl(y) is a simple polytope, and so A~ x y > for all a G Ja- 
Recall that A a G Z mxn and A/ G Z mx ( n ~ m ) stand for [A]fce<x and 
[Afe]fc£<r> respectively. Similarly, given a vector x G Z n , the vectors x CT and 
ay stand for (xk)kea and {xtjk^a respectively. The following result is from 

m P . sis]. 
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Proposition 1. Let y G M. m and let tt(y) be as in (1.1\) . and let y G 7 for 

some chamber 7. Then, 

(3.5) My)} = Yl {CMy),*)}, 

modulo the group generated by the characteristic functions of convex polyhe- 
dra which contain affine lines. 

Proof. Using notation of [HJ p. 817], define the linear mapping p : W 1 — » 
R m with p(x) = Ax, so that the polyhedra -Pa(?/) and Q(y) are identical. 
Moreover, for every basis a G B(3a,i)i v <j '■ ^ m ~~ > ^ n is the linear mapping: 

y^ [v*(y)]v = A- X y and [v (y)]j = 0, y G M m . 

Finally, for every x G M n with x > 0, p a (x) := x — v a (Ax) satisfies, 

[/v(£)]<j = -A~ l Ajij and [^(x)]^ = x^. 

Therefore, the cone [v a (y) + p a (C)] in [8] is the set of points x G M m such 
that Xtf > and x a = A^ 1 (y — A^x^); and so this cone is just [C(fl(y), a)] 
in (|3.4p . Therefore a direct application of (3.2.1) in [51 p. 818] yields (|3.5p . 

□ 

Theorem 2. Lei y € Z m , 2 G C™ mi/i \\z\\ < 1, and let y G 7 /or some 
chamber 7. Recall the set of bases B{3a->1) defined in h3.1\) . With P := 
£l(y), the rational function h defined in h 1.2(1 can be written: 

(3.6) h(y;z) = £ f[C(fl(y), a), z] = £ ^T^T ■ 

x-i * — , 5 a (y — A^u) zH 

(3.7) mtt z » Rx^z) := ^ * £ ^ . ^ 

(3.8) and z ^ R 2 {<J;z) := ]J [l - (z^-^*)^ . 

The pair {R\,R2} is well defined whenever z G C n satisfies 7^ and 
z k 7^ •2o-" 4 ' T ^ fe for every basis a G 3 a which does not contain the index k G" a. 



Proof. By a direct application of Brion's theorem to the sum (|3.5p . the 
associated rational functions f[fl(y),z] and /[C(fi(y), 0"), z] satisfy: 

(3.9) h(y,z) = m(y),z] = £ /[C(ft(y), a), z]. 

<reB(JU,7) 

Therefore, in order to show (|3.6j) , one only needs to prove that the rational 



function R ^.^ is equal to f[C(fl(y), a), z], i.e., 
(3.10) ^(^i£) = J- {z x I x G C(Q(y),a) D Z"}, 
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on the domain D a = {z £ C n | 1 > \zkZ~ A<y Ak \ for each k G" a}. Notice 
that 



1 n 



7 



II E A a <A k E ii„ . I T 1 \ rr r ' 

k&a V k eN LZ <J J V £f$n-m z a 



on D a . On the other hand, according to (|3.4p . the integer vector x G Z n lies 
inside the cone C(P, 14-) if and only if : 



x CT = A a l (y - Ajxj), 5 a (y - Ajxj) = 1 and 
xj = u + //(jU, with ti G Z™J m and v G N n ~ 



But from the definition ([32]) of a; z) and z* = zJ ' z%° = z^'^ A5/av) , 

(3.H) R j v ; a ^ = z yy y, E 5a{y ~ A ' u)z7 



— mn. — m. ~ . r— RJ r». — m. 2r/T 



which is exactly (|3.1U|) . Notice that ay = it + /vi>, and so 5 a {y — A^u) = 
5 a (y — Ajx,/) because of the definition (12. 3p of 5 a . Finally, using (13.111) in 
(|3.9p yields that (|3,6p holds whenever ||z|| < 1 and Ri(y,a;z) and R2(cr;z) 
are all well defined. □ 



Notice that i?2 is constant with respect to y, and from the definition (|2.3p 
of S a , R\ is quasiperiodic with periods and fj, a , meaning that 

^ 3 12 , #i(y + 4r<? 5 cr;^) = ^i(y,o";^)4 and 

Ri(y + fiaq,o~,z) = Ri(y,a;z) (za a q Y" 

hold for all y,q G Z m . Obviously, the more expensive part in calculating 
-R2O in (|3.8p is to compute the determinant fi a = | det -A CT |. On the other 
hand, computing Ri(-) in (|3.7|) may become quite expensive when /x CT is 
large, as one must evaluate /U™ _m terms, the cardinality of Z™~ m . However, 
as detailed below, a more careful analysis of (|3.7p yields some simplifications. 

3.1. Simplifications via group theory. From the proof of Theorem [21 
the closed forms (|3.7|) - (|3.8|) for Ri(-) and i?2( - ) are deduced from ([3.111) . 
i.e., 

Ri(y,a;z) _ ^ y ^ 5 a (y - A^ X(/ ) z x / 



v ' J ™ ./~ r Wn. — m. 



after setting ay = u + fi a v and recalling that ^(y) is a periodic function, 
i.e., 5 a (y + fi^q) = <5 CT (2/) for all y,q E Z m . However, we have not used yet 
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that S a (y + A a q) = 5 a (y) as well. For every a G Sa, consider the lattice : 

(3.13) K a := QAjZ C Z m , 

jeer 

generated by the columns Aj, j € a. The following quotient group 

(3.14) g a := Z m /A ff = Z m /0A,Z 

jeer 

= {Ec[0, a] , Ec [2, a] , . . . , Ec[/^ - 1, a] } 

is commutative, with = | det A a \ elements (or, equivalence classes) Ec[j, a], 
and so, Q a is isomorphic to a finite Cartesian product of cyclic groups Z^ , 
i.e., 

Qcr — ^"r/i X ^JJ2 X • • • X Z^ s . 

Obviously, /j, a = r/ir/2 • • • Tj s , and so, Q a is isomorphic to the cyclic group 
Z Mct whenever \i a is a prime number. Actually, Q a = {0} whenever fj, a = 1. 
Notice that the Cartesian product Z J)1 x • • • x Z^ can be seen as the integer 
space 7L S modulo the vector r/ := (7/1,772, ■ ■ ■ ,Vs)' € N s . 

Hence, for every finite commutative group Ga, there exist a positive integer 
s ff > 1, a vector r/ CT G N Sct with positive entries, and a group isomorphism, 

(3.15) ga ■ Qcr ^ V° mod r^, 

where g a {^) mod r/o- means evaluating [g a (0]k mod [r/ CT ]fe, for all indices 1 < 
k < So-. For every y £ Z m , there exists a unique equivalence class Ec^cr] 
which contains y, and so we can define the following group epimorphism, 

(3.16) h a ■ Z m -> Z Sct mod 7? a , 

y i-> My) : = 5<r(Ec[7 y ,cT]). 

On the other hand, the unit element of Q a is the equivalence class Ec[0, a] 
which contains the origin, that is, Ec[0, a] = {A a q \ q £ Z m }. 

Hence, h a (y) = if and only if there exists q S Z m such that y = A a q. 
We can then redefine the function 5 a as follows, 

^ y » == { l ±1:,°- 

One also needs the following additional notation; given any matrix B € 

Z mxt , 

(3.18) h a (B) := [h (7 (B 1 )\h a (B 2 )\---\K(B t )} G Z^ x *. 

And so, from (|3.7p . — A^u) = h a (y) — h a {A^)u mod rj a . Finally, 
using (|3.17p in (|3.7p . one obtains a simplified version of Ri(-) in the form: 



(3.19) Ri(y,a;z) 




u <E Z"- m : 



^(y) = K{Aj)u mod 7? CT 
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Next, with q G Z' m fixed, v q A a l q G Z m for some integer u q , if and only if 
v q h a (q) = mod r) a . If we set v q = fi a , then \i a A~ x q G Z m , and n a h a (q) = 

mod n a , because Q a has \i a = | det A a \ elements. Nevertheless, may not 
be the smallest positive integer with that property. So, given o~ G Ja and 
k £ a, define v k(J > 1 to be order of h a (Af : ). That is, is the smallest 
positive integer such that Vk,ah a (Ak) = mod r] a , or equivalently : 

(3.20) Vh,aK x A k G Z m . 

Obviously Vk,a < M<r- Moreover, is a multiple of z^ jCT for it is the order 
of an element in Q(j. For example, the group Z^ modulo tj — (^ has 14 
elements; and the elements b\ = (J), 62 = (?) and 63 = (*) have respective 
orders : 2, 7 and 14. Notice that, 2&i = 7&2 = 14&3 = mod rj. But, 
26 3 = 26 2 # and 76 3 = 61 ^ mod 77. 

The important observation is that o~ CT (y — f/c i(T ^-fco) = 5 a {y) for all g G Z m 
and k ^ cr, which follows from (|3.20p and ()2.3|) . Thus, following step by step 
the proof of Theorem [21 we obtain: 

Corollary 3. Let y G Z m , z G C n with \\z\\ < 1, and let y G 7 /or some 
chamber 7. Recall the set of bases B(Ia,^) defined in \3.1\) . With a G 
^(Ja,7), ^et i?i and i?2 fre as m Theorem^ Then 

r o 2n Ri{y,cr;z) _ R*(y,a;z) 

1 " J ^ 2 (y;z) i*I(y;*) ' 



(3.22) 



where : 



,0"; z 



n[> 



-A- L A k \Vk,< 



(3.23) i?t(y,a;z) := ^ » £ 



uj G E/^; 

^a(y) = h a (A^)u^ mod t? ct 



, Act Arftirf 
I Zcr 

twtt Iff := {ty G N n_m I < u k < u k><7 - 1; k <£ a}. 
One can also obtain (I3.21D by noticing that: 

Rt(y,a;z) R* 2 (a;z) ^ + P +P h 



where f3 k>a = z^z a Aa Ak , and \x„ is a multiple of v klJ . 



3.2. Simplifications via finite number of generators. Decompose Z m 
into n a := |detA<j| disjoint equivalent classes, where y,£ G Z m are equiv- 
alent if and only if S a (y — £) = 1. For every basis a € Ja, let Q a be the 
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quotient group denned in (|3.14j) . that is, 

G a = Z m /QAjZ = {Ec[0,a],...,Ec[» a -l,a}}. 

jeer 

Notice that |/,(GZ" belong to Ec[j, a] if and only if A~ l (y - £) £ Z n , and 
that Z m is equal to the disjoint union of all classes Ec[j, a]. 

Next, pick up a minimal representative element of every class, i.e., fix 

(3.24) a] G Ec[j, a] such that A~ l y > A' 1 ^, a] > 0, 

for every y G Ec[j, a] with A~ 1 y > 0. The minimal representative elements 
ct] in (|3.24p can be computed as follows: Let d G Ec[j, <r], arbitrary, and 
let d* € Z m be such that his £>entry df. is the smallest integer greater than 
or equal to the fc-entry of —A~ 1 d. The vector a] defined by d + Ard* 
satisfies (I3T24]) . 

Notice that d* + > 0. Besides, let d,y G Ec[j,cr] with A" 1 ?/ > 0. 
There exists q G Z m such that y = d + A a q. Hence q > —A~ 1 d; in addition, 
q > d* follows from the above definition of d*, and so A~ l y > d* + A~ 1 d > 
0. 

Therefore, the vector a] := d + A CT d* satisfies (I3.24p . In particular, if 
Ec[0, a] is the class which contains the origin of Z m , then £[0, a] = 0. Notice 
that for every integer vector y G Z m , there exists a unique £ [j, a] such that : 

y = eb>] + A a q, for q G Z m . 

Moreover, the extra condition A~ 1 y > holds if and only if: 

(3.25) y = £b>] + A(j q with q G N m . 

We obtain a compact form of /i(y; z) when y G Z m n 7, for some chamber 

7- 

Theorem 4. Let /i and a] 6e as in and \3.2J$ , respectively. Let y G 
Z m ri7, /or some chamber 7. Recall the set of bases B{Jai 7) defined in \3. 1)) . 
For every basis a G S(A,7) t/iere is a unique index < a] < fi a such 
that y is contained in the equivalence class Ec[j[y, a], a] defined in \3.14\ ), 
and so: 

(3.26) £ Mtfejlk^l z ^ yS _ 

CT GB(A, 7 ) 2(0", 

where \_A a ~ 1 y\ G Z m is suc/i i/tai ais k- entry is the largest integer less than 
or equal to the k-entry of A~ x y. 

Proof Recall that if y G Z m n 7 

i?l(y,cr;^) 



%;z) = V 



Next, recalling the definition (|3.1|) of ,B(Ja,7), At 1 ^ > for every basis 
a G £>(Ja,7) with y G 7. Recall that there is a unique index j[y,a] < fi a 
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such that y = £\j[y, a], a] + A a q with q G N m ; see (|3.25p and the comment 
just before. 

To obtain the vector q G N m , recall that the minimal representative ele- 
ment a], a] in (|3.24p is the sum y + A a y* where y* G 7L m is such that 
his &-entry y\ is the smallest integer greater than or equal to —A~ 1 y, for we 
only need to fix d = y in the paragraph that follows (|3.24p . In particular, 
[Aj^yJ = — y*, and a], a] = v — A„[A„~ 1 v\, which when used in (|3.7p 
and fling]) , yields, 

Rl(y,a;z) = R 1 ^[j[y,a],a},a; z) z\ A "^. 
And so (USD implies (13T261) . □ 

Theorem H] explicitly shows that it suffices to compute Ri(v,a;z) for 
finitely many values v = a], with a G £>(A, 7) and < j < n a , in order 
to calculate h(y;z) for arbitrary values y G Z m 1I7, via (|3.26p . 

In other words, in the closure 7 of a chamber 7, one only needs to consider 
finitely many fixed convex cones C(f2(£[j, cr]), a) C M ra , where cr G B(A,j) 
and < j < (A a , and compute their associated rational function (|3.26p . The 
counting function h(y;z) is then obtained as follows. 

Input: y G Z rn Hj, z G C n . 
Output p = z). 

Set p := 0. For every a G £>(A,7) : 

• Compute £[j[y,a], a] := y - A a [A- l y\ G Z m . 

• Read the value Ri(£,\j[y, cr], cr], cr; z)/R2(cr; z), and update p by: 

fli(gb[y»g].g]>g;*) j A -i y j 

^2(0"; ^) 

For the whole space Z m it suffices to consider all chambers 7 and all cones 
C(n(t\j,a]),a) C M n , where cr G B(A, 7 ) and < j < 

Finally, in view of (13,7p - (|3.8p . the above algorithm can be symbolic, i.e., 
z G C m can be treated symbolically, and p becomes a rational fraction of z. 



4. Generating function 

An appropriate tool for computing the exact value of h(y; z) in (|1.2|) is 
the formal generating function H : C m — > C, 

n 

(4.1) s -> H(s) := £ %; *) = 1^^' 

j/ez m fc=i 

where s y is defined in (12. 4p and the sum is understood as a formal power 
series, so that we need not consider conditions for convergence. This gener- 
ating function was already considered in Br ion and Vergne [8] with A = ln(s). 

Following notation of [8j p. 805], let < x G K n be a regular vector, 
i.e., no entry LA" 1 ^],,- vanishes for any basis a G Sa or index 1 < j < m. 
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Define : 
(4.2) 



1 if [A- x Ax]j > 0, 
-1 if [A^Axlj < 0. 



Next, for every basis a G I a, index j G a and vector G 
(4.3) 

Define also the vector r][a, iy] G Z n by : 



fix 



9\j,a,Uf] G Z : the smallest integer greater 
than or equal to — ej :(7 [A~ 1 Ag'U r /]j. 



(4.4) 



«j if j cr; 



1 - 0[j',<r,*y] if j G cr, Ej )(7 
The following expansion can be deduced from [8]. 



1. 



Theorem 5. Let < x G M n 6e regular and consider the vectors r][a,u^] G 
Z ra defined in \4-4\ ) f or °~ G Ja "ty G Z n ~ m . T/ie following expansion 
holds: 



(4.5) 



n 

fc=l ' 



1 



E 
E 



where 



(4.6) 



{0, 1, . . . , fj, a — 1}, /v = I det A a \ and: 
R2(cr;z) 



z 



Proof. From Brion and Vergne's identity [HJ p. 813], 



(4.7) 



n 1 r 



F(C? + P(T (C),L), 



where F(C? + p a {C),L) is the formal power series E^ ; e z added over all 
elements I in the intersection of the cone C% + p a (C) with the integer lattice 
L = Z[wi, . . . ,w n ]. Moreover, the coefficients £j )(T are defined in f|4.2j) and 
the cone C? is defined by the following formula [8j p. 805], 



(4.8) 



C% 



Finally, given the real vector space W = M>[u>i, . . . ,w n ], every p a : W — > 
W is a linear mapping defined by its action on each basis element of W, 



(4.9) 



Pa(w k ) := W k - [A^AkljWj. 

jea- 
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Hence, p (7 (wj) = for every j G a, and the cone p a {C) is given by 
(4.10) Pa (C) = 



k^cr 



jetr 



ay>0 



see [U p. 805]. Thus, every element in the intersection of the cone C? + p a (C) 
with the lattice Z[w\, . . . , w n ] must be of the form : 

(4.11) E x * Wfc + J^ejvr&wj, with x^N"~ m , 

& 6 Z m and ^ > -Sj^A-^x^j. 
On the other hand, for every basis a, define p a = | det A a \ and : 

(4.12) = + firVj, with M^Zj;;™ and ^ e N n_m . 

Moreover, as in (|4.3p . fix 9[j,a,u^] € Z to be the smallest integer greater 
than or equal to —Ej^[A~ l A^u^\j. Thus, we can rewrite (|4.1ip so that the 
intersection of the cone C? + p a {C) with the lattice Z[u?i, . . . , ui n ] must be 
of the form : 

(4.13) ^2 [ u kWk + v k PaP(wk)} + ^2 £ h° w i a ' u <Ai + 1j\ ' 

fcgcr jGcr 

with tyeZ™;"\ d^M"-" 1 and q a eN m . 

We can deduce f)4. 13|) from (|4,lip by recalling the definition (|4.9p of p a (w k ) 
and letting : 



Since F(C? + p a (C), L) is the formal power series ^ e* with summation 
over all elements I in (|4.13p . one obtains 



(4.14) 



F(C? + p a (C),L) 

E 



11 1 _ e £j,o-Wj 
jGcr 



fcger 



e ^i tT p tT {w k ) 



With 77[o-,iy] G Z n being as in using (|4.14p into (f4T7|) yields the 

expansion 

1[T-^ E £ 



- 1 

n 

3=1 



3=1 



n — - 



P<j{v>k) 



Finally, we defined w k := In (2%) +]n(s)A k for every index 1 < k < n, 
where the vectors s, z S C n have all their entries different from zero and 
ln(s) is the [1 x n] matrix defined in (|2.4p . So e Wk = z k s Ak . Moreover, 
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recalling the definition (|4.9I) of p (T (wk), the following identities hold for all 

1 < k < n, 

(4.16) Pa (w k ) = ln(z fc ) - ^ln(^)[^^fc]i- 

jeer 

Notice ^2j &a ^j[A~ 1 Aj < .]j = A k . A direct application of (|4.16p and the 
identities e Wk = z k s Ah yields (|4.5p . i.e.: 

JL i z r,[<7,u?] s A v [a,us) j_ 

1 - z k s A * ~ ^ ^ Rzi^z) Al ! _ z 



Aw ' 



fe=l 



wii/i R2(cr;z) = |l 



□ 



A direct expansion of (|4.5|) yields the following: 



Theorem 6. Xei 0<i £ R n 6e regular, and let h and ij be as in hi. 2\) and 

respectively. Let J a be the set of bases associated with A. Then for 
every pair of (y, z) E Z m x C n with \\z\\ < 1: 

(4.17) h(y;z) 



E 



' R 2 (o-,z) 

aa A) A~ 1 y>0 



V 



fU—m Z, 



A a 1 Aju 



1 ifA-^y-An^u]) eN™, 
otherwise, 



where: 




(4.18) 


< 


(4.19) 


and 



-A- L A k \Ha 



R 2 (a;z) := ]J [l - (z k z a 

k£cr 

Proof. Recall the expansion of H(s) as a formal power series : 



(4.20) 



n 

H( S ) = Hy;z)s y = Ur^Y 

y ez m k=i k 



ZkS Ak 



We also have a similar formal power series for the product : 



n i _ z .gAj n 



E ^ 



^ zfs A ^ 
Combining the latter with (14. 5j) yields: 

(4.21) = E E E 



i? 2 (o-; z) 
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Notice that (14.201) and (|4.2ip are identical. Hence, if we want to obtain 
the exact value of h(y; z) from (14.21[) . we only have to sum up all terms with 
exponent Ar][a, u<f] + A a q a equal to y. That is, recalling that each A a is 
invertible, 

(4.22) h(y-z) = E E 



Riicr: z) 



X 



1 iiA'^y-Ariia,^]) e N" 
otherwise. 



On the other hand, setting q a := A a 1 (y — Ar][a,u^]) and recalling the 
definition (fOj) . 

7?[<r,ty] CT ty A^y A~ x y «</ 

(4 23) z l\°Mz*' = - = - J- 

Ztj Ztj 

We finally prove that the vector A~ l Arj[a,u^] is bounded, so that (|4.18p 
holds, 

(4.24) < [A' 1 Ar][a,u 9/ }} j < 1 for each jeer. 
Prom (|4.3p and (|4.4p . the following identity 

(4.25) [A'JA^, u Aj - [A-MM, + { , I Z ' -1 

holds for every index j e o. 

Next, suppose that \A~ X A^u^}j = R + r where R e Z and < r < 1 are 
the respective integer and fractional parts. We can obtain 6[j, a, ty] in (|4.3p 
as follows: If £j >a = 1 then 9\j,a,u^] = —R is the smallest integer greater 
than or equal to —R — r, and so (14.25P yields 

[A~ l Ari[a,u l/ ]] j = r. 
If Sj )ff = —1 and < r < 1, then 6\j, a, Uj\ = R + 1 and : 

[A- 1 A V [a,u (/ }] j = r. 
At last, if £j >a = — 1 and r = 0, then cr, ly] = -R and so, 

[A-^r/[a,^]] . = 1. 

Therefore (14.241) holds because < r < 1. Finally, a direct application of 
(f4~23|) into (14221) yields the following version of (|4"T7j) : 



' x 



%;*) = r — ^ y 

V 7 ^ R 2 (cr;z) ^ A- X A^ 

cGjyi , A& y>0 uj&l~ m z <y 

x< 1 HA-^y-Ania,^]) e N* 
otherwise. 
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Notice that the first sum is calculated only on those bases a € J a which 
satisfy A~ l y > 0. This follows from (j4.22j) combined with (j4.24j) . □ 



Observe that (|4.17|) is different from ()3.6f) or (|3.21j) because in (|3.6p and 
(|3.2ip . the summation is over bases a in the subset B(A,y) C {JU; A~ x y > 
0}. 
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